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We report on recent results of a high statistics lattice calculation of the unpolarized and polarized struc- 
ture functions of the nucleon. 
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1 Introduction 

Our knowledge of the deep-inelastic structure functions 
of the nucleon, both for unpolarized and polarized beams 
and targets, has improved a lot in the last couple of 
years 1 . But all efforts of our experimental colleagues 
are in vain if we do not succeed in computing the distri- 
bution functions from first principles. 

This is basically a nonperturbative problem. The 
tool to solve it is lattice gauge theory. We have seen 2,3 
that this is possible with the help of dedicated computers, 
and in the next couple of years we hope to make real 
progress in this endeavor. 

The theoretical basis of the calculation is the opera- 
tor product expansion, which relates the moments of the 
structure functions to forward nucleon matrix elements 
of local operators: 
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The Wilson coefficients are computed perturbatively, and 
the reduced matrix elements v n ,a n and d n , which are 
defined by 
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are computed on the lattice. Here we have restricted our- 
selves to the leading quark operators. The structure func- 
tions are independent of the subtraction point pt. The 
moments have the parton model interpretation: 



v (J) = {x n-i ) (J) ; a ^ = 2 Au,a 
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where x is the fraction of the nucleoli momentum car- 
ried by the quarks and Am, Ac? denote the quark spin 
contribution to the spin of the nucleon. The structure 
function consists of two contributions: a^/- 1 has twist 



two, whereas 
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has twist three. The latter contribution 



has no parton model interpretation. 

2 Highlights of the Calculation 

The (euclidean) lattice operators are obtained from (3) 
by the substitution 
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where is the lattice covariant derivative 



Choosing 

suitable interpolating fields B, B for the nucleon, the de- 
sired nucleon matrix elements can be derived from ratios 
of three-point to two-point functions, 
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for < t < t. 

The calculations in this talk are done on 16 3 x 32 lat- 
tices at [3 = 6.0. We use Wilson fermions and work in the 
quenched approximation, where the effect of dynamical 
quark loops has been neglected. To be able to extrapo- 
late our results to the chiral limit, and to make contact to 
the predictions of the quark model, we do the calculations 
at three different hopping parameters, k = 0.155,0.153 
and 0.1515, which correspond to physical quark masses 
of roughly m q = 70, 130 and 190 MeV, respectively. So 
far we have analyzed of the order of 1000, 800 and 400 
independent gauge field configurations, respectively, at 
the three values of k. 

The chiral limit is obtained by extrapolating m 2 lin- 
early to zero in 1/k. For the critical hopping parameter 
we find k c = 0.15693(4). At our smallest quark mass, 
taking the nucleon mass as the scale, the inverse lattice 
spacing is a -1 k, 1.4 GeV, which corresponds to a dis- 
tance of 0.14 fm and a spatial box size of 2.2 fm. In 
the chiral limit the inverse lattice spacing varies between 
1.8 and 2.2 GeV, depending on how we extrapolate the 
nucleon mass to the chiral limit. 

The operators must be chosen such that they belong 
to an irreducible representation of the hypercubic group 
if (4). We consider the operators 
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They require at most one spatial component of the nu- 
cleoli's momentum to be non-zero 2 . 

The lattice operators must be brought into a contin- 
uum renormalization scheme, e.g. the momentum sub- 
traction scheme or the MS scheme. Schematically this 
is written 

O(fi) = Zo(afi)0(a). (7) 

We have computed the renormalization constants 
Zq perturbatively to one loop order 2 and non- 
perturbatively 4 following a recent suggestion in the lit- 
erature 5 . The non-perturbative calculation uses O(100) 
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Figure 1: The renormalization constants Za lQ ,Zy 2 and Zy 3 as a 
function of jj? (in lattice units) for n = 0.1530. The solid (dotted 
and dashed) lines are the results of standard (tadpole improved) 
perturbation theory as discussed in Ref. 3. 
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Figure 2: The moments (x), (x 2 ) and (x 3 ) as a function of 1/k. 

gauge-fixed lattice gauge field configurations. In Fig. 1 
we compare both results for some of our operators. We 
find that for larger values of fi 2 perturbative and non- 
perturbative renormalization constants agree within 20% 
and better. At smaller values of fi 2 we cannot really 
expect to find agreement with perturbation theory yet 
at our value of [3 because here the internal momenta are 
no longer large compared to any hadronic scale. In the 
following we have taken the perturbative values for the 
renormalization constants. The numbers we will quote 
are for fi = a -1 . 

3 Unpolarized Structure Functions 

Let us first discuss our results for the unpolarized struc- 
ture functions. The moments (x" -1 ) are shown in Fig. 2 



for our three quark masses. The lines are linear fits to the 
data, and the solid circles are the results of the extrapo- 
lation to the chiral limit. The lowest moment is found to 
be practically independent of the quark mass, while with 
increasing n the moments show a growing increase with 
the quark mass. 

One is tempted to compare the quenched quark dis- 
tribution functions with the phenomenological valence 
quark distribution functions. But this is not quite cor- 
rect because we are dealing with probabilities. Consider 
(x), and let the subscript V, S and G denote the valence 
quark, sea quark and gluon contribution. Momentum 
conservation then demands that (x)^ + (x)y^ + (x)s + 
(x)g = 1. In the quenched approximation (x)s = 0, so 
that valence quark and gluon contributions should sum 
up to one. If, for example, the gluon contribution is ap- 
proximately equal in both cases, we would thus expect to 
find that the quenched quark result is noticeably larger 
than the valence quark contribution. A better quantity 
to compare is the difference oft/- and rf-quark distribution 
functions. In the chiral limit we obtain - {x)^ = 

0.26(4). Martin, Roberts and Stirling 6 (fit £>_) find 0.18, 
while the CTEQ collaboration 7 (fit CTEQ3M) obtains 
0.22. Our result is one standard deviation away from the 
CTEQ result. For the higher moments we find agree- 
ment with the phenomenological numbers within the er- 
ror bars. 

4 Polarized Structure Functions 

We shall now turn to the disscussion of the polarized 
structure functions. Let us first concentrate on Am and 
Ac?. In Fig. 3 we show our results together with a linear 
fit to the data. The solid circles indicate the extrapola- 
tion to the chiral limit. For heavy quark masses we find 
Am r; 1 and Ac? as —1/4; m good agreement with the 
quark model 8 . In the chiral limit we obtain 

Am- Ac? = g A = 1.07(9). (8) 

This difference may be expected to be insensitive to the 
effects of dynamical quarks. The experimental value for 
the axial vector coupling constant is g A = 1.26. Our 
result and the experimental value differ by two standard 
deviations. 

If we combine this result with the perturbatively 
known Wilson coefficients we obtain for the Bjorken in- 
tegral 

/' dx(g{(x, Q 2 ) - g1(x, Q 2 )) = 0.174(25) (9) 
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at Q 2 = a~ 2 as 3 — 5 GeV 2 . Because of uncertainties in 
the mass extrapolations to the chiral limit as mentioned 
before, we cannot state the scale more precisely at the 



moment. Our result (9) is in agreement with the exper- 
imental value 9 of 0.163(17). For the second moment of 
gi we find 

/ dxx 2 (g p 1 (x,Q 2 )- g1(x,Q 2 )) = 0.0162(54). (10) 
Jo 

This number is also consistent with experiment. 

Let us finally discuss the structure function g^ . An 
interesting question is how large the contribution of the 
twist-three operator is. In Fig. 4 we have shown c?2 for 
our three values of the quark mass together with a linear 
fit. In the limit of heavy quark masses c?2 goes to zero 
as one would have expected. In the chiral limit, on the 
other hand, c?2 is compatible in magnitude to 02 (which 
we have no space to show here 4 ). Altogether we obtain 

/ dxx 2 (g p 2 (x,Q 2 )- g? 2 (x,Q 2 )) = -0.0257(60). (11) 
Jo 

Without the twist-three contribution the number would 
have been -0.0109(35). 

5 Conclusions and Outlook 

We have seen that an accurate calculation of the struc- 
ture functions of the nucleon from first principles is fea- 
sible. Our results so far are encouraging. If they do not 
exactly reproduce the experimental values this may have 
various reasons. 

So far we have only considered the leading quark 
operators and neglected the contributions of gluonic op- 
erators. Gluonic operators are much harder to get at. A 
preliminary study of matrix elements involving gluonic 
operators shows however that it is possible to compute 
the gluon distribution functions accurately as well. We 
hope to present first results in the near future. 

For a quantitative analysis at low Q 2 it will also be 
important to include the contribution of higher twist op- 
erators. As we have seen in the case of c?2, higher twist 
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Figure 3: The quark spin contributions to the spin of the nucleon 
Au and Ad as a function of 
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Figure 4: The twist-three matrix element d,2 as a function of 

contributions may be large, contrary to one's expecta- 
tions. 

Other sources of uncertainties are finite cut-off and 
volume effects as well as the lack of dynamical quark 
loops. Certain tests of finite cut-off effects have shown 
that they are presumably small 2 . We know how these 
uncertainties can be systematically reduced, and we will 
do so in the future. 
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